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Abstract 

Convection-diffusion problem are the base for continuum mechan- 
ics. The main features of these problems are associated with an indef- 
inite operator the problem. In this work we construct unconditionally 
stable scheme for non-stationary convection-diffusion equations, which 
are based on use of new variables. Also, we consider these equations 
in the form of convection-diffusion-reaction and construct uncondi- 
tionally stable schemes when explicit-implicit approximations are used 
with splitting of the reaction operator. 



1 Introduction 

Convection-diffusion equation are basic in the mathematical modelling of the 
problems of continuum mechanics. The main features of these problems are 
connected with the nonseldadjoint property of elliptic operator and dom- 
ination of convective transport. When considering compressible media, an 
operator of convection-diffusion problem is indefinite. In this case, given pro- 
cess can be nondissipative, i.e. norm of the homogeneous problem solutions 
does not decrease with time. This behavior of the norm solutions need to 
pass on the discrete level in choosing of approximations in time. 

In the numerical solution of non-stationary problems for convection-diffusion 
equations the most widely used two- and three-layer scheme. Investigation 
of the stability and convergence of approximate solutions can be performed 
using the general theory Samarskii A. A. of stability (correctness) of operator- 
difference schemes 011]. Must be kept in mind that for convection-diffusion 
problems direct application of the general stability criteria can be difficult 
due to non- self adjoint operators. Note also that, in view of indefinite operator 
of problem we need to oriented ^-stable (g > 1) operator-difference schemes. 
In the solution non-stationary problems of long periods of time preference 
should be given asymptotically stable schemes [5]. For these schemes ensures 



1 



the correct behavior of the solutions with the release of the fundamental 
solutions for large time and damping of others. 

In this paper, we construct unconditionally stable scheme for the approxi- 
mate solution of non- stationary convection-diffusion problems. Such schemes 
can be applied to other problems with an indefinite operator. The study con- 
ducted by the example of a model two-dimensional boundary- value problem 
in a rectangle. Used the simplest approximation of the operators of diffusive 
and convective transfer on a uniform rectangular grid. Constructed uncondi- 
tionally ^-stable difference scheme based on the introduction of new variables 
and the explicit-implicit approximation. 



2 The convection-diffusion problem 

We consider the Neumann problem in a rectangle for the non-stationary 
convection-diffusion equation. For simplicity, assume that the coefficient of 
diffusion transport is a constant (independent of time, but depends on the 
point of the computational domain). The coefficient of convective transport 
is natural to consider the variables both in space and time. 
In the rectangle 

Q = {x \ x = (xi, x 2 ) , < x a < l a , a — 1, 2}. 

We consider the non- stationary convection-diffusion equation with the con- 
vective transport in divergent form, 

=/(*>*)> 0<t<T, 

a=l ^ ' 

in the standard assumptions ki < k (x) < k 2 , k\ > 0, T > . This equation 
is supplemented by Neumann boundary conditions 

du(x,t) xe dSl, 0<t<T. (2) 

on 

For the unique solvability of the nonstationary problem the initial condition 
is given 

u(x,0) = u°(x), xett. (3) 
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On the set of functions u (aj, t), which satisfy the boundary conditions (j2J), 
non-stationary convection-diffusion problem written in the form of differential- 
operator equation 
du 

— + Au = f(t), A = C(t)+V, 0<t<T. (4) 
The diffusion operator T> is defined by 

^ -til: (*<»>!:) 

and convection operator C 

2 d 

a=l a 

Cauchy problem is considered for the evolution equation (Pll) : 

u (0) = u°. (5) 
For convection operator we have the following representation 



1 1 2 / 

C = C + -divv£, C u = -^2(v a (x,t) 

a=l ^ 



du d , , . . 
+ tj — {v a (x,t)u 



where £ — the identity operator and Cq -the operator of convective transport 
in a symmetric form. 

For arbitrary functions u(x) n w(x), we define the Hilbert space % = 
L2 (fl) with inner product and norm 

(u,w) = / u(x)w(x)dx, \\u\\ = (u, u) 1 ^ 2 . 
Jq 

Diffusion operator D on the set of functions satisfying (J2J), is self-adjoint 
and positive define 

V = V* > 0. (6) 

The operator of convective transport is considered under the assumption that 
the normal component of the medium velocity v = (ui,^) on the boundary 
is zero: 

v n (x) = vn = 0, x G <9f2, (7) 

where n — outward normal to dVL. In H the convection operators have the 
following properties: 

Co = -C *. (8) 
Also useful upper estimates for convective transport operator C: 

\(Cu,u)\<5\\u\\ 2 , S = ^\\dwv\\ c{n) . (9) 
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3 The differential-difference problem 

For an approximate solution of the non-stationary convection-diffusion prob- 
lem we use a uniform grid in the area Q: 

u = {x | x = (xi, x 2 ) , x a = + ^ h a , 

V = 0,l,...,jV a) (N a + l)h a = l a , a = 1,2}. 

We define the Hilbert space H = L 2 (oo) for grid functions, where the inner 
product and norm are defined as follows: 

G/> w) = ^y(x)w (x) hih 2 , \\y\\ = (y, y) 1/2 . 

For the difference operator of the diffusion transfer D is used additive 
representation 

2 

D = J2 Dia \ a = 1,2, sew, (10) 

Q=l 

here D^ a \ a = 1, 2 is associated with the corresponding differential operator 
in one direction. 

The difference operator of diffusion transport ffTQ]) in is self-adjoint and 
positive definite [3] 

D = D*>0. (11) 

The convective terms are approximated with second-order, using the cen- 
tral difference derivatives and shifted grids to specify the velocity compo- 
nents. For the difference operator of convective transport are also using 
additive representation 

2 

Cy = Y,C {a) . (12) 

a=l 

For the case of sufficiently smooth velocity components and solutions of 
the differential problem, for example, we can assume, 

b a (x) = v a (x), x G ft, < x a < l a , 

t J — ■ x Q — - x — / q j ot — X . _ . 

Difference operator of convective transport in symmetric form have the 
following basic property: 

C * = -C . (13) 
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We also |S] have the grid analogue of inequality 

\{Cy, y )\<5\\ y \\ 2 

with a constant 



(14) 



r 1 

o = — max 
2 



(xi + 0.5/ii, x 2 ) - (a* - 0.5/ii, x 2 ) 



+ 



6 (2) (xi, x 2 + 0.5/i 2 ) - b (2) (xi, x 2 - 0.5/i 2 ) 



h 2 

Therefore, from the equation (jlj) we arrive at the differential-operator 
equation 

^ + Ay = <p {t) , A = A(t) = C + D, 0<£<T, (15) 

on the set of grid functions y(t) e H with the initial condition 

2/(0) =2/°. (16) 

Difference convection and diffusion operators in the differential-difference 
problem inherit the basic properties of differential operators. 



4 Unconditionally stable schemes 

For simplicity, we restrict ourselves to a uniform grid in time 

u T = u T U{T} = {t n = nr, n = 0, 1, N, tN = T}. 

For an approximate solution of (ITS]) . ( fl6|) commonly used two-layer weighted 
scheme, which have a following restrictions on the time step r < r = 

To construct the unconditionally stable schemes for the solution of the 
differential problem (115"]) . (Tl6|) with ^4 > — SE, 5 > we define a new function 

UK 

y = exp(St)w. (17) 

Substitution of (TTT1) in (TToT) . (Tl6|) with homogeneous right-hand side gives the 
following problem for the w: 

^ + Aw = 0, A = A + SE, 0<t<T, (18) 

w(0) = y°. (19) 
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Under this transformation, operator A is a nonnegative (A > 0). 

To solve the problem (ITS]) . (lH?j) we use a two-layer weighted difference 
scheme, which is unconditionally stable for standard restrictions a > 0.5. 
We write the scheme for the grid function y n 

ex V {-5r)y n+1 - y n + ^ + ^ ^ ae ^_ 5 ^ y n+i + (1 _ a ) y *) = , (20) 

y° = u°, t n E 0J T , ■ (21) 

In contrast to the non-standard schemes considered in [6], the positive effect 
is achieved not only through the use of a new approximation of the time, but 
also by correcting the problem operator. 

Theorem 1. The difference scheme / fl?0j) . I[21\) with a > 0.5 unconditionally 
g-stable in H with 

g = exp(5r), (22) 
with the a priori estimate for solutions 

\\y n+1 \\ <Q\\y n l (23) 

Proof. We rewrite the scheme (1201) . f[2"Tj) in form 

^W-^ +i/+ l = 0i t n eUr) (24) 

T 

where 



p n+l 



aexp(-Sr)y n+1 +(l - a) y n = r (a - r n+l + l - (exp(-dr)y n+1 - y n ) , 
r n+i = exp(-5T)y n+1 - y n 

T 

Multiplying the scalar equation (T2~^|) by p n+1 , we obtain the equality 

(r n+ \r n+1 )+A(p n+1 ,p n+1 ) 



+ ^ ({exp(-5r)y n+1 ,exp(-6r)y n+1 ) - (y n ,y n )) = 

From this equation, under the condition a > 0.5 and A > 0, yields the 
estimate of stability ([23}, ([22]). □ 
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Equation ([T]) can be written in the form of convection-diffusion-reaction 
equation with the convective terms in the symmetric form 

du 

— + C u + Du + Ru = f(x,t), xett, t>0. (25) 

\J v 

where ^ 

Ru = r(x,t)y, r(x, t) = — div v. 

For the reaction operator we have the estimate 

R = R*, -58 <R< 5E. (26) 

To construct the unconditionally stable scheme without the assumption 
of nonnegativity operator of problem we will use the explicit-implicit ap- 
proximation for the equation (125]) . The problem is generated by the reaction 
operator therefore we split it into two: 

R = R+ + R_, R+ = R* + , f?_ = fT, 0<R+<5E, -5E<R_<0. 

(27) 

When using the two-layer explicit-implicit schemes, we can only count 
on first-order accuracy in time. Therefore it is natural oriented to purely 
implicit approximation of the basic terms of the operator and define following 
difference scheme 



71+1 

'' ' ~ +(C n + D + R n + )y n+l + R n _y n = 0, n = 0,l,...,N- 1. (2* 



T 

Theorem 2. Explicit-implicit difference scheme |I2J), ( tUP unconditionally 
g-stable in H with 

q=1 + St (29) 
for the numerical solution we have the estimate: 

h n+1 \\ < Qh n \\, n = 0,l,...,N-l. (30) 

It is important to note that, in contrast to the ordinary weighted scheme, 
stability is obtained without restrictions on the time step. The transition to 
a new time layer associated with the solution of the grid problem 

(E + ar(A + 5E))y n+1 = X n (31) 



for scheme ( 120]) and 

{E + r{C + D + R + ))y n+1 = r n (32) 

for scheme (}2"H|) . 

Equation (13"T|) and ( I32p is a system of linear algebraic equations with 
a positive definite nonself adjoint matrix. For solutions that can apply the 
standard iterative methods. 
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